Behavior of a frustrated quantum spin chain with bond dimerization 
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We clarified behavior of the excitation gap in a frustrated 5 = 1/2 quantum spin chain with bond 
dimerization by using the numerical diagonalization of finite systems and a variational approach. 
The model interpolates between the independent dimer model and the 5=1 spin chain by changing 
a strength of the dimerization. The energy gap is minimum at the fully- frustrated point, where a 
localized kink and a freely mobile anti-kink govern the low-lying excitations. Away from the point, 
a kink and an antikink form a bound state by an effective triangular potential between them. The 
consequential gap enhancement and the localization length of the bound state is obtained exactly 
in the continuous limit. The gap enhancement obeys a power law with exponent 2/3. The method 
and the obtained results are common to other frustrated double spin-chain systems, such as the 
one-dimensional J1-J2 model, or the frustrated ladder model. 
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I. INTRODUCTION 



Recently, the interest in the field of the low- 
dimensional quantum systems is concentrated on systems 
with a spin gap, especially, in connection with the high- 
Tc cuprates upon doping.EI The ground state of such a 
system realizes the spin disordered singlet state, favored 
by a quantum fluctuation and/or frustration. The ex- 
citation energy gap opens above the ground state. 
typical models, we can consider the spin ladder model J3'EI 
the boHd-alternation model,H and the Majumdar- Ghosh 
model.Qu „ 

Syntheses of various new compoundsO also accelerate 
investigations on this field both theoretically and experi- 
mentally. For example, magnetic susceptibility measure- 
ment on KCUCI3Q indicates a spin gap behavior, and the 
experimental data are considered to be explained theoret^ 
ically by the double spin-chain model with frustration.u 
Estimation of the magnitude of the gap through the sus- 
ceptibility data is |already established on the assumption 
of the dimer gap.El However, it is not sufRcient to the 
thorough understandings of the system. We must inves- 
tigate the origin of the gap and the physical picture of 
the excited states. This is a main purpose of this paper. 

The A chain, the subject model in the present paper, is 
a participant of double spin-chain systems, but has rather 
special geometry of the interaction bonds. The trian- 
gles are aligned in one direction. If we connect the spins 
located at the top of the triangles with the interaction 
bonds, it becomes the railroad-trestle model. Therefore, 
it may seem that the interesting properties are peculiar 
to this special system. 

The early investigations on the A chain were directed 
in seacch-jfor a system with the singlet dimer ground 
state.ll3~EJ Interpretation of the excited staljes by a kink 
and an antikink was done at this stageJlj'EJ The exis- 
tence of the finite pjergy gap above the ground state was 
rigorously proven.Ej Analogy of the model to the kagome 



antiferromagnet has been pointed out:0 there are macro- 
scopic local continuous degeneracy in the ground state in 
the classical limit, low-lying excitation spectrum is con- 
sequently dispersionless, and there exists an additional 
peak of the specific heat at low temperatures. The origin 
of the dispersionless mode and thej-deuble peakj-of the 
specific heat were recently clarificd.ESEj Sen et aJi3 also 
pointed out that the possible relevance of the model to 
the newly synthesized compound, YCUO2.5. 

In this paper, we find that the A chain possesses com- 
mon features with other spin gap systems besides its 
unique properties stated above. Concretely, we mean the 
"common feature" by the way how the system reduces 
frustration by the bond dimerization. This is an essence 
of this paper. 

The ground state is generally unstable against a small 
perturbation that relaxes strong frustration. It must 
have a strong infiuence on the low-lying excited states as 
well, and consequently to the low-temperature behavior 
of various physical quantities. We should take this effect 
into account when we analyze the experimental data. We 
consider the bond dimerization as a perturbation, since 
its energy stabilization is the strongest one, and it can 
be realizjed by a lattice distortion as in the spin-Peierls 
systems. t3 

From a theoretical point of view, the frustrated double 
spin-chain systems with the bond dimerization are very 
attractive. They have the dimer state and the S = \ Hal- 
dane state in both extremes of the strength of the dimer- 
ization parameter. In the midst is the fully-frustrated 
point. As was shown in the bond-alternation model,u 
there always exist the dimerj-ardeillS and the string order 
of den Nijs and Rommelse.llj The dimer state continu- 
ously changes to the S = 1 Haldane state with the in- 
crease of strength of the ferromagnetic interaction bonds. 
Therefore, there is always a finite energy gap. Then ques- 
tions arise. How do we identify these two phases? How 
the differences manifest in the observable physical quan- 
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titles? In this paper, we show by adopting the A chain 
for an example that the whole phase space can be di- 
vided into three regions with respect to the nature of the 
first excited state relevant to the energy gap, and that 
the differences affect the low-temperature dependences 
of physical observables such as the specific heat and the 
susceptibility. 

Section || describes the model and summarizes the gen- 
eral remarks concerning the symmetric A chain without 
the dimerization, which will be the starting point of the 
discussion in Sec 
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In the first part of Sec 
do the same variational analysis as was done in the sym- 
metric A chain.E3E3 This is valid only for positively (anti- 
ferromagnetically) small dimerization. Then, continuous 
limit of the effective Hamiltonian is derived and its exact 
The remaining part of Sec 
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solution is obtained, 
devoted to the case when the dimerization is negatively 
(ferromagneticallvL small. We use the non-local unitary 
transformations Ej to represent the ground state and 
the excited state. This transformation is equivalent to 
the Kennedy- Tasaki transformation of the S = 1 system, 
and is its adaptation to the double spin-chain systems. 
It is known to be powerful when the ground state is ei- 
ther in the Haldane state or in the state with strong 
dimer correlation. Almost equivalent results to the pos- 
itive dimerization case are obtained near the symmetric 
point. We show how the system converges to the Haldane 
state and the dimer state in Sec. Section shows 
the differences of the observable quantities between two 
phases. We also propose a quantity to judge the phase 
in experiments. Section VI is devoted to summary and 
discussion. 



II. THE MODEL AND GENERAL REMARKS 



other limit of A = — oo, the system becomes equivalent 
to the S = I spin chain. Therefore, the present model 
has the independent dimer ground state and the Haldane 
ground state in its extremes. 

The understand™^ of the symmetric A chain are al- 
ready established .lla'ta Here, we briefly summarize the re- 
sults. The ground state is the perfect singlet dimer state 
with two-folded degeneracy in the case of the periodic 
boundary conditions. The low-lying excited states ap- 
proximately consist of (A^ — I) singlet dimers and two free 
spins. The excitations are governed by sets of two free 
spins named a "kink" and an "antikink" . A kink stays 
localized and works as a delimiter to moving antikinks. 
Dispersionless aspects of the excitations originate in a 
localized kink. An antikink is considered as a free par- 
ticle moving between kinks with the effective mass. An 
antikink is supposed to be one free spin only within the 
first approximation, since it is not an eigenstate of the 
local Hamiltonian. It spreads out to an extent in reality. 
A detailed structure of an antikink is not revealed yet, 
but it merely renormalizes the effective mass. Within 
the first approximation, the mass m = 1. The second 
approximation that an antikink is distributed among 5 
spin sites gives the mass m — f .158,c3 and the numeri- 
cal diagonalization data after the N ^ oo extrapolation 
shows the mass m = 1.2I.L3 The excitation gap is well 
expressed by a summation of a kinetic energy of an an- 
tikink and the creation energy of a pair of a kink and an 
antikink, eq = 0.215, as 
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We consider a system described by the following 
Hamiltonian. 



N 



= ^ XcTn ■ Tn + Tn ■ CT „+i + (T „ • CT „+i 

n=l 



(1) 



Here, A^ is the number of the triangles in the system, A is 
a parameter denoting the dimerization, and \(t\ = |t| = 
1/2. Figure ^ shows the depicted lattice. 
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FIG. 1. Shape of the dimerized A-chain. Bold lines indi- 
cate the A bonds. 



where k is the wave number of an antikink alone. The 
low-temperature peak of the specific heat can be repro- 
duced by the rSchottky specific heat with the gaiP-| ex- 
pressed above.113 Susceptibility was also calculated.ll3 



III. IN THE VICINITY OF A 



At the point of A = 1, the system is the symmetric A 
chain and is fully frustrated. In the limit of A = -l-oo, the 
system reduces to the independent dimer model. In the 



We numerically diagonalized the above Hamiltonian up 
to the systems with 28 spins (A^ = 14) under the periodic 
boundary conditions. 
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where denotes a singlet dimer state connecting the 
ith and the jth site, namely = (tiij — Ii1j)/V2 for 
Ti (ii) denoting an up (down) spin located at the ith site. 
Here, we numbered the site so that cr„ is the {2n — l)th, 
and Tn is the (2n)th as is shown in Fig. |^. The t2i above 
is an antikink. The wave function of a kink located at 
the leftmost edge is known as, 



V'kink = [Tl (T2i3 - i2T3)+ T2 (Tlia - ilT3)]/V6. (4) 



FIG. 2. The A dependence of the energy of the lowest state 
in the subspace denoted by k (wave number), and R (spin 
reversal symmetry) in the system with A'^ = 14 (28 spins). 

The A dependence of the energy is shown in Fig. |^. 
The ground state energy is highest at the fully-frustrated 
point (A = 1), which leads to the instability due to the 
lattice dimerization. 

For A > 1, the ground state is always the pure dimer 
state which consists of the singlet dimers on the every 
A-bond. Thus the total ground state energy is exactly 
— 0.75AA^. We refer to this state as the left-dimer state 
hereafter in this paper. In this region, the spin corre- 
lation length vanishes. The excited states also have no 
fc-dependence. At A = 1, the ground states are two- 
fold degenerate under the periodic boundary conditions. 
They are the left-dimer state and the right-dimer state. 
The first-order transition occurs at this point. For A < 1, 
the correlation length gradually increases with decreas- 
ing A. The right-dimer state smoothly reaches the S = 1 
Haldane state in the limit of A = — cxd without any phase 
transition. Hereafter, we simply call the region A > 1 as 
the dimer phase, and the region A < 1 as the Haldane 
phase. The degeneracy of the excited states in the dimer 
phase is lifted in the Haldane and a state with fc = tt 
becomes the first excitation. 



A. A > 1: the dimer phase 



This state is an eigenstate of the edge triangle Hamilto- 
nian, and therefore does not move. Its energy eigenvalue 
is A/4 — 1. Thus a kink contributes to the excitation 
energy by (A — 1). The singlet dimers, {[2n, 2n + 1]}, ex- 
isting between a kink and an antikink, are not eigenstates 
of local triangle Hamiltonians. They also contribute to 
the excitation energy. 

A variational basis is not orthogonal to each other and 
satisfy the following relations. 



(V'jIV'j) 
{i'^\n\yJJ) 



E^ + iX- 1) 



+ - (A - l)min(i,j) 



-0,;. 



(5) 
(6) 



Here, Sij is the Kronecker delta, and min(i, j) = J if i < j. 



Our task is to find a function ^'var 
imizes the energy expectation. 



C'iT/'i that min- 



Ya.TE 



(*var|W|*v 



(7) 



The results -obtained in this subsection has been re- 
ported briefly.cfl We discuss the method and the results 
in detail in this paper. 

The excited states in this phase may have similar 
properties to those at the symmetric point, because the 
ground state is the same. Therefore, we do the-same 
variational analysis as was successful at A = I MM The 
boundary conditions are set open, though it does not in- 
fiuence the final results in the thermodynamic limit. We 
define a variational basis ipi so that an antikink is located 
at the ith triangle. 

= V'kink ® [4, 5] • • • [2i -2,2i- 1] 



If we diagonalize the denominator and rewrite the nu- 
merator with its eigenfunction, this variational problem 
is transformed into a simple eigenvalue problem. Then, 
the solution of Eq. (0) can be obtained by the numerical 
diagonalization for a finite system size. We show the re- 
suh of iV = 200 in Fig. | with A = 1.00, 1.001 and 1.01. 
In the case of the symmetric A chain (A = 1.00), the wave 
function is the sine function indicating a free motion of 
an antikink. An antikink is drastically attracted to a kink 
as A increases. The wave function exhibits an antikink 
localized. We clarify the wave function analytically by 
using the continuous limit. 
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FIG. 3. Variational wave function of an antikink in the 
lowest excited state for A — 1.00, 1.001, and 1.01. Size of the 
system A'^ — 200. n stands for the site number of the triangles. 



We rewrite the above equation by the new basis \4>k) = 
exp[iA:n] |i/)„), since the denominator of Eq. (|^) is di- 
agonalized by Fourier transformation in the large N limit. 
Note that k is the wave number of an antikink alone, and 
does not correspond to the total wave number. Then the 
basis relations become 



h\n\(j)i) 



k,l 



5 — 4 cos k 

= [E^ + (A - 1): 

+ ^(A-l)Afc,i, 



(8) 



(9) 



(i5g + (A-l) + ^ + y 14,- 
13 1 

^iV4^ '{k-i + isy 



(12) 



Apart from the constant terms, this Hamiltonian is 
equivalent to the following in the real space represen- 
tation. 



— — -(\ 



l)xexp[— (5x]. 



(13) 



Here, x is the distance between a kink and an antikink; 
m is the effective mass of an antikink and is set m — 1 
in the first approximation. We take the limit ^ — > at 
this stage. The first term is the kinetic energy of an an- 
tikink, and the second term is the triangular potential 
attracted by a localized kink. We rescale x hy X — 6x 
with 9 = (3m(A — l)/2)~^/'^. Then the eigenvalue equa- 
tion Tic^ — Eq'^ becomes 



dX2 



(14) 



where E' = Ec x 2m6^. Its solution is known as the 
Airy function, = Ai{X — E') ■wpili the first eigenvalue 
E' ~ (37r/2 X 0.7587)2/3 ~ 2. 338.^1 Accordingly, we can 
obtain the energy eigenvalue Ec, the average distance 
between a kink and an antikink (x) , and the localization 
length of the wave function ^. 



with 



Ak,i = ^exp[z(A:n - ^m)] min(m,n) ('(/'„ 1^/;™)- (10) 



The diagonal elements of Akj dominate the off-diagonal 
elements. Then, Ak k is written as 



1 



6 — 13 cos fc -I- 8 cos 2k 



N 245 
4 



- 87 cos k - 
4 — 5 cos k 



30 cos 2k + 4: cos 3fc 



^ - 10 cos /c 
N + 1 3 



2 cos2fc 



2 5 - 4 cos A: 



(11) 



We find in this equation that Ak_k diverges in the limit 
N — !■ oo. In order to avoid this divergence, we introduce a 
cut-off factor iS to the momentum, namely k k-\-i6/2. 
This should have no effect on physical results by taking 
the limit (5^0 after N ^ oo. We rewrite the matrix 
element and pick only up the leading term of the off- 
diagonal part and the terms up to the fc^ in the diagonal 
part. Then the continuous limit, N ^ oo and fc — > 0, of 
the Hamiltonian T-lk,i is given by. 



En — 



E' 
2"^ 



= 1.532m-i/3(A- 1)^ 



x) = -E'e = lM2m-^'^{\ - l)-i 
3 

^ ~ 5 X 6* = 4.368m"^/^(A - l)"^. 



(15) 

(16) 
(17) 



The localization length,^, is obtained from a rough es- 
timation of the localization length of the Airy function 
^ 5. The estimation of ^ is quite consistent with the 
wave function shown in Fig. |: ,^ - 40 for A = 1.001, 
and ^ ~ 20 for A = 1.01. It should be noticed that 
jAk^k/ {(l>k\(t>k) of Eq. (^) and (^ becomes equivalent to 
Eq, if we replace N by the ^ above. 

Now, the total gap behavior can be obtained from the 
kink contribution, (A — 1), and the antikink contribu- 
tion, Ec- The gap enhancement defined by Acap — 
£^gap(A) -£;gap(l) is 

Acap = (A - 1) + 1.5319 X m-^'^{\ 



N2/3 



(18) 



The gap increases in a power law with its exponent 2/3. 
Figure ^ shows the behavior of the gap enhancement, 
Acap, compared with the numerical results of the peri- 
odic system of iV = 14. 
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FIG. 4. Log-log plot of the gap enhancement for the exact 
result in the continuous limit with the mass m — 1.21, and 
the numerical results of the periodic system with A*' = 14. 

Wa-tised a value of the mass, m = 1.21, estimated at 
A = llla for better comparison. For A — 1 > 0.01, the nu- 
merical data agree with the analytical estimation. In this 
region, the localization length is within the finite system 
size, ^ < 14. 

The relevant excitation in this region is governed by a 
competition between the kinetic energy and the potential 
energy of an antikink. When the dimerization is small, 
an antikink gains energy by the kinetic motion. As the 
dimerization becomes large, an antikink is bound by a 
kink, and finally they collapse to a local triplet. The 
exponent 2/3 is a general outcome of this competition. 



B. A < 1: the Haldane phase 

In this region, the ground state cannot be known triv- 
ially, although we can expect the right-dimer state con- 
tinuously changes |-tn— tie S = 1 Haldane state in the 
limit of A ^ — ooJl3St3 We make use of the non-local 
unitary INLIJ) transformation for the double spin-chain 
systemsjH3c3 the second-order perturbation, and the nu- 
merical diagonalization to clarify the ground state first. 
Then we proceed to investigate the excited state. 

We transform the Hamiltonian Ti. with U defined in 
Appendix ^. 



N 

E 

ri=l 



U-^nU = > : Acr„ • x„ 



- 4(« + fT,^T^X+l<+l- (19) 

According to this transformation, string order parame- 
ters are transformed into local parameters such as (cr^), 
(cr^), (cr^r^), as will be found in Eqs. and (|2j). It 
makes possible that we employ a single-site approxima- 
tion like the molecular-field one. Thus, we consider the 
following variational trial function for the ground state. 



I*0> = n = 11(^1^") + VT^I^n)) (20) 

n— 1 n— 1 

|5„) = (|T,i>-U,T»/V2 (21) 

|T„) - a| T, T) + /3(| T. i) + I i. T))/V2 + 7| i, i) (22) 

I t, T)'s are the state of |cr^, r^). 6, a, 7 are the varia- 
tional parameter and satisfy the normalization condition, 
+ + 7^ = 1. These parameters are supposed to be 
invariant of n. The analysis is variational in this sense. A 
state with 6 = is the left-dimer state on the A bonds, a 
state with h — V3/2 is the right-dimer state, and a state 
with b = 1 corresponds to the pure VBS state in this rep- 
resentation. The energy expectation value is calculated 
as 



(*o|W|*o) 



A b 



2\2 



+ 2/32(a2+^2) 



36Vl - ^'¥(«' - 7') 



N. 



(23) 



We can easily find this minimum value and the varia- 
tional parameter set by using the Lagrange multiplier. 
The energy value eg is 



2&V3(l-&') _ 



(24) 



with four possible choices of the parameters {a, (3,^) as 

(25) 



(a,/?,7) = (±V2/3,VlA0), 
(0,-7173, ±7273), 

and b determined through 



b{Ab^ - 3) 



0. 



(26) 



(27) 



The four- fold degeneracy in the choice of a, /L^x-tpr re- 
sponds to the degeneracy of the edge states.ElfLJCJ A 
state with 6 = is the singlet dimer ground state for 
A > 1. The other one, Eq. (p6|), corresponds to the 
ground state for A < 1. If we solve Eq. ( p6| ) up to the 
second order of (1 — A), 
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1 



1- A 



:(1-A)^ 
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Then the energy expectation per triangle is 



(28) 



(29) 
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This agrees with the result due to the second-order per- 
turbation of (A — 1)<T„ • T„. Details of the calculation 
are given in Appendix ^ The diagonalization result of a 
finite system with 14 triangles shown in Fig. || is fitted 
by the least-square method to 



FIG. 5. The string order parameter and the dimer order pa- 
rameter are plotted against A/(l — A) for A < 0, and against 
A/(l -|- A) for A > 0. Bold lines are the analytic estimate 
obtained by the non-local unitary transformation. Thin lines 
are numerical diagonalization results of the lattice with 12 
triangles. 



E/N = --- 0.1886(1 - \f - 0.003(1 - A). (30) 

Consistency between the analytic estimation (note 
3/16 = 0.1875) and the numerical result is excellent. 

We can also estimate the steing order parameter of 
den Nijs and Renwielse, Ostr,li^ and the dimer order 
parameter, Odim-"l23 The definitions and the expectation 
values are 



dim 



lim 

\m — 'n\—f 



-4([/ V^exp 



k—m-\-l 



hm 4(a>,t>=4«.)«>=4(0' 

\ra — 7i\—*oo 

lim -{U-'S:,,exp[zn ^ Sl]S^U) 

1 7n — n\—>-oo ^"^^ 

k—m-\-l 

hm (5^5^ 



(31) 



z \2 



ml 



(32) 



At A = 1, Odim = 1 as & = \/3/2. In the limit of 
A ^ — oo, it converges to the VBS value Odim = 4/9 as 
6^1. On the other hand, the string order parameter, 
Ostr takes a value of Ostr = l/4atA = l and converges 
to the same value of Odim — Ostr = 4/9. It should be 
noted that the convergence of 6 ^ 1 is very fast that it 
takes a value h = 0.94 even at A = 0.25. Therefore, the 
Haldane state can be realized even if all the interaction 
bonds are antiferromagnetic. 

Figure |^ shows the dimer and the string order param- 
eter of the ground state in the whole phase space. 
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Numerical data are of the periodic system with N = 
12. Quantitative agreement is excellent near A = 1. The 
string order increases as A decreases from 1. Therefore, 
we call the region A < 1 the Haldane phase. It should 
be noted that thepHumerical data converges to the value 
of 5* = 1 (~ 0.37)c3 in the limit of A ^ — oo, not to the 
value of the pure VBS state. The contradiction becomes 
distinct at A ~ —1, where the string order parameter 
takes the maximum value. This is because we used the 
variation of the single-site approximation as in Eq. (|2^) , 
besides the correlation length is rather large in the Hal- 
dane state. 

A low-energy excitation is intrinsically of domain-wall 
type in one dimension. Two perfect singlet dimer states, 
|n(\/3/2)) and |«(0)), are degenerate at A = 1. For A < 1, 
the latter state becomes the excited state and the former 
state becomes the ground state by changing the value of 
h. Therefore, we consider the following wave function as 
a trial for the excited state. 



N 



\n—l n—i-\-l 



(33) 



Because the state |n.(0)) becomes a higher excited state 
for A < 1, the validity of this trial wave function is re- 
stricted to the very vicinity of A = 1. The domain wall 
located at the ith triangle is essentially an antikink be- 
fore the NLU transformation is done. Thus, the analysis 
is for a kink-antikink excitation. The basis relations of 
ipi are 



mi:,) = (Vl- 62)1-^1 (34) 
{^^\n,\^J) - (-Eg + EiiN- min(i, j)) + Sj^,,,) 

m^^j), (35) 

with = eoA^, El = -0.75A - eo, and E2 = Xib"^ - 
0.75) — Eq. These relations are equivalent to the ones in 
the previous subsection, Eq. (|^) and (^. Therefore, the 
problems are solved in the same way and we obtain 



Ec = 1.856m"i/3£:f , 
{x) = 1.22,lm-^/^E~^ 
^ ~ 3.969m-^^^Ep 



(36) 
(37) 
(38) 



for the excitation of /c = 0. However, we remark that 
the lowest excitation is the state with fc = tt in this re- 
gion. This state is considered as that both a kink and 
an antikink are mobile with the total momentum tt. The 
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analysis on the state of A: = tt will be done in the next 
subsection. Three estimates, Eqs. (|3^), (|37|), and (|3^ ) 
above, however, have quite good consistency with the 
second excited state with fc = 0. 



IV. AWAY FROM A = 1 

Away from the fully- frustrated point, A = 1, the exci- 
tation is no longer of a kink-antikink type. They already 
collapse to a local triplet state. Therefore, we must con- 
sider another type of a domain wall for a trial wave func- 
tion. 

In the region of A < 1, the ground state has four- fold 
degeneracy associated with the edge states. Within the 
scheme of our variational analysis, it appears in the four 
possible choices of the variational parameter (a,/3,7) of 
Eq. (^). The most natural candidate for the domain 
wall is the one between any two of the four-fold-degener- 
ate ground states. In fact, Fath and J. Solyomcj showed 
that the lowest excitation in the AKLT modeS is of this 
type. The trial wave function we consider is then 



i?cx(fc) always takes minimum at fc = tt. The energy gap 
Ecx{t^) converges to 1/9 in the VBS limit, A —oo. 

We can also discuss excited states of A > 1 with this 
variation scheme, since the domain wall of Eq. ( p9| ) can 
express the local triplet domain wall, if we change the 
notation of the pair, (T„ and t„. Unless, the variational 
solution for the ground state only gives the state with 
6 = 0, since the exact ground state is the left-dimer 
state on A bonds. Then, the variation with 6 = gives 
nothing at all. It should be noticed that a basic recipe 
of the present variational method is that a spin pair, 
Sn — Cn + Tn, should be chosen so that it does not 
take a singlet dimer state. Therefore, we shift cr spins 
by one site as cr„ — s- (t„+i. With this new definition 
of cr„ and t„, the left-dimer ground state is represented 
by 6 = \/3/2, or in other words, a = 0. Expressions of 
El and E2 become different from Eqs. (42) and (Q), and 
are Ei = —A/2 — 1/4 and E2 = A, respectively. The 
excitation is calculated dispersionless as 



EaAk) = A. 



(45) 



N 



\^i)=Y.Q^..=J2cAl[\n{b)) n |n'(6)) 



\7i—l n—i-\-l 



(39) 

where |n(6)) and \n'{b)) only differ the choice of the set 
(a,/3,7). We use the set {a,P,j) = {y/2/3, ./T/3,0) for 
|n(6)), and (a,/3,7) = (-^273, v/T/S, 0) for |n'(6)). 6 is 
common and determined by Eq. ( pq ) . The basis relations 
are calculated in the same way as 



(^.IV',) = (1-/ 



[E, + {\i-j\ 



^(-a)l^-^"l, 
l)Ei+S,^,{Ei+E2)] 



(40) 



with 



Ei^--b^ 1 6+V3(l-62) + 

3 I (5+^3(1^)3 



E2 = -63(6-1- v/3(l- 62)). 



(41) 



, (42) 



(43) 



In the thermodynamic limit, the above matrices are diag- 
onalized by the Fourier transformation. An expectation 
value of the energy gap above the ground state is esti- 
mated as 



E,x{k) = 



\n\ 



-eJi- 



2a {l + a'^)coak + 2a 
1 — a2 1 + 2a cos k + a'^ 



1 + 2a cos fc + 



E2 



1 + 2a cos k + a 



I -a' 



(44) 



This result is nothing but the local triplet excitation, 
where one singlet dimer is replaced by a triplet in the 
ground state. This becomes the exact solution in the 
limit A — > 00. 
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FIG. 6. A dependence of the energy gap obtained by 
the numerical diagonalization (circles), the variation of a 
kink-antikink type (broken line), and the variation of a lo- 
cal triplet type (solid line). The second excitation gap with 
fc = (triangles) is also plotted for A < 1. a:: axis is A/(l — A) 
for A < 0, and A/(l + A) for A > 0. 

Figure || shows the A dependence of the energy gap 
obtained by the variation and the numerical diagonaliza- 
tion of the periodic system with 24 spins. The numerical 
results are depicted by symbols. The variational results 
of a kink-antikink type excitation, Eq. ( |l5| ) and (|3^), are 
shown by a broken line; those of a local triplet excitation, 
Eq. (^) and (^5|), are shown by a solid line. We also plot 
the second excitation gap for A < 1 for the comparison 
with Eq. (|6|). 
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In the dimer phase, A > 1, consistency between the 
numerical and the variation is exceUent. Estimations by 
the kink-antikink variation is better in the vicinity of 
A = 1, while those of the local triplet becomes better as 
A increases. This crossover occurs at A ^ 3, where the 
average distance between a kink and an antikink (x) of 
Eq. dig ) becomes unity. On the other hand, the con- 
sistency in the region, A < 1, only remains in a qual- 
itative level. Within the single-site approximation em- 
ployed in this paper, cr spins and r spins are equiva- 
lent to each other. For example, (n(6)|cr„ • (Tn^i\n{b)) is 
equal to {n{b)\Tn ■ T„+i|n(6)). Therefore, the approxi- 
mation should be better for the system with a symmetry 
of exchanging cr ^ t. Since the present model does not 
possess this symmetry, the estimation is not good. It 
should be excellent in the.symmetric systems such as the 
Majumdar- Ghosh modelE We must go beyond the single- 
site approximation to improve the estimates in the Hal- 
dane phase of the A chain. The numerical data converges 
to the value, ~ 0.24 consistent with the 5 = 1 system; a 
half the gap of a typical diagoaalization result of S" = 1 
system with 12 spins, 0.4842.c3 Only the number of in- 
teraction bonds between Sn and Sn+i determines the 
strength of the effective interaction, since cr„ and t„ be- 
comes symmetric in the S = I limit. The present model 
has only two interaction bonds, which corresponds to a 
half the effective interaction in the 5 = 1 system. 



V. OBSERVABLES 



We calculated the magnetic susceptibility and the spe- 
cific heat of finite system {N = 7) with periodic bound- 
ary conditions in order to see the difference of observable 
physical quantities between the dimer phase and the Hal- 
dane phase. In other words, we try to determine the 
phase by these observables. 



0.20 



>? 0.15 




Figure ^ shows the uniform magnetic susceptibility for 
A = —5.0, 1.00 and 3.0. These parameters correspond to 
the Haldane phase, fully-frustrated point, and the dimer 
phase, respectively. We rescale the temperature by each 
energy gap in order to see the qualitative differences, 
namely E {E — Eg)/Ef,ap- The peak width and the 
height among three are quite different. Data of the dimer 
phase have a sharp peak, while it becomes broad in the 
Haldane phase. This reflects the structure and the den- 
sity of the excited states, i.e., many continuum states of 
multiple-magnon excitations in the Haldane phase bring 
about a broad peak. Contrary, excitations in the dimer 
phase are considered to be discrete local triplet excita- 
tions, which generates a rather narrow peak. Full width 
at half maximum (FWHM) of the peak for each data in 
the unit of T/Gap is: - 2 for A = 3.0, - 4 for A = 1.0, 
and ~ 6 for A = —5.0. The FWHM is almost three times 
as wide in the Haldane phase as in the dimer phase. We 
speculate that this value might be a judge to determine 
the phase. Recently, we have also calculated the sus- 
ceptibility of the frustrated ladder model, and the Ji- 
J2 model with bond dimerization (J1-J2-J3 model), and 
found that the FWHM take almost the same |jj"alue as the 
present case. Details are reported elsewhercu 
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FIG. 7. Uniform susceptibility per spin calculated for 
A = 3.0, 1.0 and —5.0 by the exact diagonalization of the finite 
systems with the periodic boundary conditions. The size of 
the system is A'^ = 7. Each A represents the dimer phase, the 
fully- frustrated point, and the Haldane phase, respectively. 



FIG. 8. The specific heat per spin calculated for 
A = 3.0, 1.0 and —5.0 by the exact diagonalization of the 
finite systems with the periodic boundary conditions. The 
size of the system is A*' = 7. 

Figure || shows the specific heat. Qualitative ten- 
dency is same as the susceptibility, except for the data 
of A 1.00 showing the double peak structurp^ which is 
characteristic of the kink-antikink excitation.E3 Data of 
the dimer phase is explained by the Schottky type spe- 
cific heat. 



VI. SUMMARY AND DISCUSSION 



We have investigated the A chain with the bond dimer- 
ization. The model interpolates the independent dimer 
model and the 5=1 spin chain model by changing A. 
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As for the ground state, the first order transition between 
the dimer phase and the Haldane phase occurs at A = 1. 
To be strict, the ground state for A > 1 is the perfect 
singlet dimer state, while that for A < 1 continuously 
changes from the dimer state at A = 1 to the S = 1 
Haldane state at A = —oo. 

As for the excited state, we can distinguish the inter- 
mediate region, in the vicinity of A = 1, from the dimer 
region and the Haldane region. This region is charac- 
terized by strong frustration. The energy gap is very 
small caused by the unstable ground state. The ground 
state has a strong dimer correlation. The excited state 
is governed by two free spin called a kink and an an- 
tikink. They become a local triplet in the dimer and the 
Haldane region. We can conclude that the intermediate 
region serves as a buffer between the other two regions. 

We clarified the wave function of the bound state of a 
kink and an antikink, and how they collapse to a local 
triplet. The essential point is that a kink and an antikink 
is bound by a triangular potential and the competition 
between this attractive potential and the kinetic energy 
causes the gap enhancement with exponent 2/3. The 
triangular potential is an outcome of unfavored singlet 
dimers between a kink and an antikink. This mecha- 
nism generally occurs when the frustration is relaxed by 
bond dimerization. For example, Chitra et al investi- 
gated the 1-dimensional J1-J2 model with the bond al- 
ternation, and found that the-sap behaves with (5^/^ with 
the alternation parameter SSB This can be explained in 
the same manner as in the present model. 

We have also calculated the uniform susceptibility and 
the specific heat to see how the phases are characterized 
by the physical observables. We clarified the qualita- 
tive difference between the Haldane phase and the dimer 
phase. The susceptibility data shows a broad peak in 
the Haldane phase, accompanied by continuum multiple 
magnon excitations. Therefore the FWHM of the peak 
may serve to judge which phase the data belong to. 
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N 



c/ - n 



Un = P++Pn exp[z^^,^], 

n-l 



P± = 



^71 



1 

2 

cr „ 



1 ± exp 



k=l 



(Al) 
(A2) 
(A3) 
(A4) 



where P+ (-P,7) is the projection operator onto states 
with the even (odd) number of Sf — ±1 for i < n — 1. 
The spin operators, a^, a^, a^, are transformed as. 



C/cr^f7 = exp 
Ua^U = exp 
Ua^U = exp 



N 

m— n+1 



fe=i 



exp 



N 

- E 

rn— n+l 



QX 



k=l 



(A5) 
(A6) 
(A7) 



In the above derivation, the following relations are uti- 
lized: 

u-^ = U, 

(exp[i7rS'"])^ = 1, 

xp± ^ X 

exp[*^^f] = - exp[z7r5^]a;: (a ^ /3), 
exp[i7rS'^] exp[i7rS'^] exp[i7rS'^] = 1 (for cychc a,/3, 7), 

which are easily obtained by using exp[i7rtTj^] = 2icr". 
Transformations for t„ are obtained by replacing a by 
r. By using the above relations, exchange interactions 
becomes. 



UcTn ■ tJJ 



' n+ 



' n+ 



lU 
lU 



a ' nj 

X X 

^x X 



A X Z ^Z Z 
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APPENDIX A: NON-LOCAL UNITARY 
TRANSFORMATION 

In this appendix, we summarize the non-local unitary 
transformation for the double 5* = 1/2 spin-chain sys- 
tems. The transformation is defined by U in the follow- 
ing. 



APPENDIX B: SECOND-ORDER 
PERTURBATION OF THE GROUND STATE 
ENERGY FOR A < 1 

In this appendix, we show that the second-order per- 
turbation of the ground state energy can be done in this 
case without knowing an excited state. We divide the 
Hamiltonian Ti into the unperturbed part Tio and the 
perturbation part Tii. H ~ Ho + Ti-i, where 

N 

Hq — E ■ Tn + ■ ""n+l + Tn ' CTn+l, (Bl) 

= (A- 1) E (B2) 

n=l 
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Here, note that the location of «•„ and t„ is different 
from Eq. (19) and Fig. |l| so that the pair S„ = o-„ + t„ 
takes the singlet dimer state for the unperturbed Hamil- 
tonian. The unperturbed ground state \ipo) is a direct 
product of singlet dimer states \Sn) on the spin pair cr„ 
and T„, i.e., 



N 



m = n 1'^") 



(B3) 



Let us first examine the first order perturbation Ei: 

N-l 

E, = (^"iTiil^") = (A - 1) ^ (B4) 



ra=l 



where 



n n 



AT 



(B5) 



|T,j) denotes a singlet state formed by two triplet pairs 
of Sn and Sn+i, i-C, 

IT,?) = (|0„)|0„+i) - |l„)|(-l)„+i> - |(-l)„)|l„+i))/x/3. 

(B6) 



|0„), |1„) and |(— l)n) denote triplet states of a pair 
Sn = fn + Tn with their eigenvalues of = (t^+ t^. 
It follows from Eqs. (^) and (psl) that 



(B7) 



since the singlet state \ S„) is orthogonal to the triplet 
state. It is found in Eq. (B5) that the operation Tn-Cn+i 
on ji/i") transforms the two singlets \Sn) and |S'n+i) into 
triplets leaving the other singlets unchanged. 

The second-order perturbation E2 is calculated as 



E2 = -{^°\ni 



1 



Ti-Q — Eq 
(A -1)2 J2 K 



1 



(A-l)^ 
Eo 



n.m— 1 

N-l 00 



^0 ^ £-0 



n.m— 



lfc=0 ^ ^ 



To proceed further, we divide Tio into a diagonal part 
and an off-diagonal part when it operates to \ipn). An 
off-diagonal part is 



K 



OD 



(c ri+l + T n+l) • C 



n+2- 



(B9) 



A diagonal part is the rest of the Hamiltonian, Tig — 
The eigenvalue of this diagonal part is Ef' = Eq + 1 = 
—3N/A + 1. On the other hand, the operation of to 



\ipi)n generates a state with three triplets located at the 
nth, the {n + l)th, and the (n -|- 2)th triangles, and these 
triplets form a singlet state. Namely, 

iV'^)=HO°iV'i)=^ n (BIO) 

where the three-triplets state 

IT^) ^ _L[(|i„^2)|(-l)„+i) - |(-l)„+2)|l„+i)) X |0„) 

+ (|l„+i)|(-l)„) - |(-l)„+i)|l„)) X |0„+2) 

+ (|l„)|(-l)„+2) - |(-l)„)|l„+2)) X |0„+i)] (Bll) 



forms a singlet state. Now we get 

no\^l,} = E?\^i) + \4,l). 

Similarly, we find 

no\^l) = E^\^Pl) + \^pf,), 



(B12) 



(B13) 



where E^ denotes diagonal energy and the is a new 
state with four triplets forming a singlet state located at 
the triangle site from n to n -I- 3. In general, we have 



fe+i 



{non^i) = iEfrm + Y.Ci\^Pl), (B14) 



1=2 

where {ipn) contains ^-fl triplets at triangle sites n ^ n+l 
forming a singlet state; and C'l is a constant. Thus we 
obtain 

{i^l\K\^n)^Km{E?)''{i^lHl.) 

^Sn,m{E?)'^. (B15) 
16 

Finally, we get the energy correction of the second order 

Eo 



Eo = 



3 (A -1)2 ^.^^^y. fEl_ 
16 Eq I p„ 



n=l fe=0 
^2 



16 



(A-ir(iv-i)^_^^ 



= -Y^(A-if(iv-i). 



(B16) 



This agrees with the variational result given by Eq. ( |29| ) 
in the text. 
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